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The response of notch filters to sudden excitations is analyzed. Unit 
step and stepped trigonometric inputs are considered for the class of 
filters derived from low-pass networks by a frequency transformation. 
It is possible in some cases to approximate the transient solutions in terms 
of Laguerre functions and deduce general properties of notch filters from 
these solutions. The use of phasing sections to modify the transient response 
is also examined. It is shown that this method can be used to effectively 
reduce the overshoot in the response to a stepped trigonometric excitation. 

I. INTRODUCTION 

In order to accurately measure noise levels in compandored com- 
munication systems, it is necessary to set the compandor charac- 
teristics at approximately the values associated with signal trans- 
mission. This is done by applying a so-called "holding tone" which is 
subsequently removed by a notch filter incorporated into the measuring 
set. This investigation originated in connection with the design of such 
a notch filter for an impulse noise counter. 1 The filter has to meet 
both frequency and time domain requirements. The frequency response 
requirements could be readily met with existing filter design pro- 
cedures. However, the time domain characteristics of notch filters 
needed investigation with regard to the suitability for the present 
application. The time domain requirement is that the filter when 
combined with a C-message weighting filter 2 and excited with a stepped 
trigonometric time function at the notch frequency should, in the 
transient state, have only a specified overshoot level. This requirement 
is imposed by the necessity to distinguish in the measuring set between 
sudden gain and phase variations and impulse noise. 

To examine the transient response of notch filters, a class of such 
filters derived from low-pass niters by a frequency transformation is 
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considered. The transient response to a unit step function, and to a 
stepped trigonometric function with the notch frequency, is expressed 
in terms of the low-pass impulse response. It is shown that the low- 
pass and notch filter response are for both excitations related by a 
Hankel transform. Some general properties of the transient response 
are deduced by considering notch niters for which the response can 
be obtained approximately in terms of generalized Laguerre functions. 

This investigation shows that notch filters would distort narrow 
time pulses of duration less than one-half the notch frequency period. 
The amount of distortion is related to the notch depth. The response 
of a notch filter to a stepped trigonometric function can be kept to a 
low level only after a certain time duration, which depends on the 
filter parameters. The transient response of a notch filter followed by 
a low-pass filter may still assume large values at relatively short times 
from the beginning of the response. A method of decreasing the 
transient response at such times by the use of phasing sections is also 
presented. The use of phasing sections is of particular importance 
where it is necessary to modify the transient response without affecting 
the frequency response. 

Although this work is primarily concerned with notch filters, the 
methods used may also be applied to determine the transient response 
of high-pass and bandpass filters, when derived from low-pass filters 
by a frequency transformation. 

II. TRANSIENT RESPONSE OF NOTCH FILTERS TO A UNIT STEP FUNCTION 

A class of notch filters derived from low-pass filters by a frequency 
transformation 3 is considered. The transformation corresponds to re- 
placing the inductances and capacitances in the low-pass filters with 
parallel and series resonant circuits, respectively. Let Tl(s) be the 
low-pass transfer function in the complex frequency domain s. The 
transformation is given by 

where w = 2ir/ , f = notch frequency, and /3 is a constant. (For 
low-pass filters normalized such that for s = ju the 3-dB bandwidth 
is at w = 1, /3 is the 3-dB circular bandwidth of the notch filter.) The 
transfer function of the notch filter in the complex frequency domain 
z, Tn{z), is related to the low-pass transfer function Tl(s) by 

M*)-n(^). (2) 
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To investigate the transient response of notch filters, that response 
is related to the impulse response of the low-pass filter. Such a rela- 
tionship is obtained by first expressing the transfer function of the 
low-pass filter in terms of the Laplace transform of the impulse re- 
sponse, jhit), 

T L (s) = f <r«f L (t)dt; (3) 

Jo 
then from (2) the transfer function of the notch filter by 

F,» =/ o "«Xp- [y£j] Um. (4) 

The time response, fusit), of the notch filter to a unit step function 
can now be obtained by inversion of the Laplace transformation, 
through integration in the complex plane over the contour T, as follows : 

where y is a positive constant. The relationship between step response 
of the notch filter and the impulse response of the low-pass filter is 
obtained by the substitution of (4) into (5) yielding 

'" = h L t /." exp " ( m. ) h(u)du dz - (6) 

In Appendix A it is shown that (6) can be expressed as follows: 
f N s(t) = \f" h(u)du]'l(t) 

~ f\ T AM J^ Ji(2Jp™)du\ Jo(2u, ^xt - x*)dx, (7) 

where 1(0 is the unit step function and J„(y) is a Bessel function of 
order n. 

In (7) the first term is the undistorted unit step response and the 
second term expresses the distortion introduced by the notch filter. 

The first term can be simplified by observing that from (3) 

J~ h(t)dt = T L (fl). (8) 

Without loss of generality, T L (s) can be normalized to be equal to 
unity at zero frequency. The expression in the brackets of the second 
term contains a Hankel transform, 4 of order one, hence (7) can be 
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written 

Us(t) = 1(0-2 r*i(2Vjte) J^/ (2co V:ri-x 2 )dx, (9) 

where ^1 is the Hankel transform. 

The integration in (9) can be performed approximately by using 
the method of stationary phase. This is a good approximation when 
the Hankel transform is a slowly varying function in comparison to 
the Bessel function. It is also shown in Appendix A that with the 
stationary phase approximation (9) can be approximated by 

f NS (t) = {i.o - ^ | o °° mu) y j^fj 1 (^d)du^nt). (io) 

An example considered subsequently shows that (10) is indeed a good 
approximation for /3 « w , i.e., for narrowband notch niters. 

III. TRANSIENT RESPONSE OF NOTCH FILTERS TO A UNIT STEPPED SINE 
FUNCTION 

The transient response to this function is obtained in a manner 
similar to the response to a unit step function. However, a more 
general response function, fwmit), is considered, 

T ( gg "\ 

2m+l r 1 L I 2 , 2 ) 

*■» - -tr /,*" (Sl&X *- (11) 

providing for the possibility of a low-pass notch combination. The 
time response of a notch filter to a stepped sine function with the 
notch frequency is obtained as a special case, by setting m = 0. 
Proceeding in a manner similar to Appendix A, it can be readily 
shown that (11) can be expressed as follows: 

f Nm (t) = u» £ I | o M f L (u) [ <J>i{i ~ x) l""j.<2VR»)«fa] 

■J m [2u <x(t - x)~]dx. (12) 

Equation (12) also contains a Hankel transform but of order m. 
Similarly as before, (12) can be approximated by using the method 
of stationary phase. It can also be readily shown, from the results in 
Appendix A, that with this approximation (12) is given by 

f Nm (t) « sin L t - ™ \ j" f L (u) ( 2 |)"' ! ; n (*)rf«-l (t). (13) 
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It is shown below that, for tn = 0, (13) is a good approximation 
for « oj„. 

IV. APPLICATION TO PARTICULAR TYPES OF NOTCH FILTERS 
4.1 Introduction 

To gain some insight into the behavior of the transient response, a 
low-pass transfer function is chosen for which the integrals (10) and 
(13) can be evaluated. As mentioned before, these integrals are Hankel 
transforms. A class of functions for which Hankel transforms can be 
evaluated in closed form are the generalized Laguerre functions. 8 - 6 In 
fact, for these functions the Hankel transforms are self-reciprocal. 5 
These functions form orthogonal sets so that, in principle, any passive 
filter impulse response can be expanded in terms of these functions. 
However, the impulse response of some low-pass filters may be ap- 
proximated closely by a single Laguerre function. Such a response can 
be used to obtain an estimate of the notch filter response. 

The low-pass transfer functions considered are of the following form : 

ll{S) ~ ( S /ai + D" +m+1 U ^ 

with m and n integers, n, m ^ 0, and «i > 0. 

The transfer function (14) is normalized such that 7^(0) = 1.0. This 
function can be considered as consisting of (m + 1) cascaded low-pass 
filter sections and, for a 2 = -«i, of n phasing sections. For physical 
readability of the n cascaded sections, it is necessary and sufficient 7 8 
that |a!2 1 ^ a\. 

The impulse response corresponding to the transfer function (14) is 9 

hn) = jrpsn* (|;)W*-'«a<*. - om (is) 

where L„ (x) are generalized Laguerre polynomials given by 10 

« w = i(-i>>(^)f-;- (i6) 

These polynomials are oscillatory for positive arguments and mono- 
tonically increasing for negative arguments. 

4.2 Unit step response 

The integral (10) with f L (u) given by (15) can be reduced to a 
tabulated integral for the special case n = 0. The applicable integral 
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is of the form 11 



Jo 



e-*V» + * +1 J„(2xVr)da: = ^ e- Y Y l '^U(Y). (17) 



For this special case the low-pass transfer function is 

1 



Tl(s) = 



(s + l) m+1 ' 



(18) 



where without loss of generality ai is set equal to one. 

The transient response of the notch filter is obtained by using (10), 
(15), and (17) and is approximately 



Mt) 



1.0 — £-«-<«*>' 



Lj,(|«)sin«o<ll(0. (19) 



It is of interest to determine the accuracy of (19) for different 
values of m. This is done in Appendix B for m = 0, 1, and 2 where 
(19) is compared with the exact solutions. It is shown that the accuracy 
of (19) is of order /3/« when expressions multiplied by both sin wt 
and cos ut are considered. However, expressions multiplied by sin ut 
only are of order (/3/w ) 2 , so that near the maximum amplitudes of 
the second term in (19), expected near cos u t tt 0, the approximation 
can be said to be of order (/3/co ) 2 . Equation (19) is therefore a very 
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Fig. 1 — Laguerre functions, e~ x L\n{x). 
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Fig. 2 — Values of a m (x). 

good approximation for 0/u> o « 1, i.e., for notch filters in which the 
3-dB bandwidth is much smaller than the notch frequency. 

In Fig. 1 the function e- x V m {x) is shown for m = 0, 1, 2. It follows 
from this figure, in conjunction with (19), that the distortion of the 
unit step function will be particularly pronounced in the time vicinity 
T = ir/2a) . In Fig. 2 the values of the second term in (19) are plotted 
as a function of 0T/2 for m = 0, 1, 2. This graph gives an estimate of 
the distortion of the unit step function, and is of particular significance 
when considering the transient response of notch filters to pulses of 
short time duration. This investigation shows that pulses of duration 
less than one-half notch frequency period will be considerably dis- 
torted by notch filters. 

To obtain a numerical estimate of the distortion, a notch filter with 
the following requirements is considered: (i) Notch frequency— 1010 
Hz. (it) Notch depth at least -30 dB in the frequency range 995-1025 
Hz. A notch filter with these requirements is derived from the low-pass 
transfer function (18) for m = 0, 1, 2, 3. 

It readily follows from the transformation (1) that /3 can be deter- 
mined from the following equations. 

/3 = o A 2*(f - h) (l + £) , (20) 

where f = notch frequency, /i is the lower frequency where a notch 
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depth A[dB] is required, and, from (18), 



"' " V exp -(io(^+l) lnl0 )- L (21) 

The 3-dB bandwidth B of the notch, also obtained from (1), is 

B = 7T~ , (22) 

where <o B is given by (21) with A = —3.0. 

Table I lists the filter parameters, including the maximum value, 
a m , of the second term in (19). 

It is evident from Table I that the 3-dB bandwidth and the distor- 
tion term a m decrease as the number of sections increases. The differ- 
ence is particularly pronounced between m = and m — 1. A further 
decrease in the distortion term can be obtained by decreasing /3 ; how- 
ever, this also reduces the depth of the notch. 

The listed values of a m seem to be representative of what is obtain- 
able with notch filters of specified notch depth and 3-dB bandwidth. 
For example, computation of a two-section notch filter derived from 
a Tschebyscheff low-pass filter with 0.5-dB ripple gave a value for 
a m of 0.18. The notch depth of this filter was also —30 dB and the 
3-dB bandwidth 170 Hz. The lower value obtained with this filter can 
be attributed to the more oscillatory behavior 12 of the low-pass impulse 
response. Equation (10) suggests such an interpretation. 

4.3 Transient response to a stepped sine function 

The low-pass transfer function (14) and the corresponding impulse 
response (15) are again considered. With that impulse response the 
integral (13) is given by 

f Nm (t) = —j- J - TWl ^J sin ^ - T j 

■ f K e- aiu u m ' 2 L?[( ai - <x 2 )u]J m (-l2J3ut)du. (23) 

Jo 

The integral (23) can be evaluated in closed form, again yielding 
Laguerre functions 6 

0t (a 2 — ai) 1 



**»-*ttWi(t) v ~'W 



2 «ia2 J 

.*,(«,( -2?). (24) 

290 THE BELL SYSTEM TECHNICAL JOURNAL. FEBRUARY 1974 





Table 


1— Notch filter parameters 




m 


"it 


COB 


0[kHz] 


B[Hz] 


dm 



1 
2 
3 


31.61 
5.53 
3.00 
2.15 


1 

0.64 
0.51 
0.43 


6.00 
1.05 
0.57 
0.41 


955 
260 

178 
149 


0.527* 
0.273 
0.234 
0.227 



* Computed from the exact expression, and occurs at a time, t - 0.12 ms. 



It is noted that with the condition for physical readability of the 
n cascaded sections in (14), |« 2 | ^ a h that the argument of the 
Laguerre functions is always positive, and hence the functions are 
oscillatory. 

From (11), (14), and (1), (24) corresponds to the inverse of the 
following Laplace transform : 



T(z) = 



~U+« a JL 



(? + o> 2 + 0/aiz) m 

Z 2 + «j 

z 2 + co 2 + fi/aiz 



] 



The terms in the brackets in (25) can be interpreted as transforms 
of (i) a stepped sine function, (ii) an w-section low-pass filter, (Hi) 
el notch filter section, (iv) n phasing sections for a 2 = -«i or n addi- 
tional notch sections for a 2 -* °° • 

The transient response of a notch filter followed by a low-pass 
filter is of interest. However, (25) is restricted to a particular low-pass 
filter with a high-frequency cutoff in the vicinity of the notch fre- 
quency, and will not be considered further. 

Setting m = 0, (24) simplifies to 



/*.«) = <-<"*->' n[£k( i -%)]**** 



(26) 



The special case a 2 -> <*> , treated previously for the unit step re- 
sponse, can be obtained from (26) and is in agreement with the result 
obtained by performing the integration directly by using (17). 

F or a2 _»oo and m = 0, 1, 2, a comparison was made between the 
exact solutions and the approximate solution (26). The comparison 
showed the same accuracies as for the unit step response. 

The time response due to a stepped cosine excitation can be ob- 
tained by differentiating (26) and dividing by w.. It readily follows 
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that, within the accuracy of (26), the same expression is obtained but 
with sin o> t replaced by cos o) t. 

It is noted that (26) gives the correct value for t = 0. This value 
can be obtained by using the initial value theorem of Laplace 
transforms. 13 

lim sF(s) = lim f(t). (27) 

Graphs of e~ x L%(x) are shown in Fig. 3 for n = 0, 1, 2. These graphs 
give the envelope of the transient response (26) for a 2 — » °° • The 
effect of finite values of a 2 can be deduced from the graphs. For 
example, for a 2 negative the arguments of the Laguerre functions 
increase reaching maximum values of (pt)/ai for ai = — a 2 . Therefore, 
for the same ftt the spacing between the zeros would decrease and 
the maximum values increase. For positive ai the opposite would be 
the case. 
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Fig. 3 — Laguerre functions, (T x L Q n (x). 
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It has been shown above that the transient response of notch filters 
can be obtained from the Hankel transform of the low-pass impulse 
response. This property can be used to deduce the qualitative charac- 
teristics of notch filters derived from conventional low-pass filters of 
the Bessel, Butterworth, and Tschebyscheff type. The impulse re- 
sponse of conventional filters, with transfer function polynomials of 
the same order and approximately the same bandwidth, has similarities 
to the impulse response considered. Therefore, the graphs shown in 
Fig. 3 are also representative of the transient response of notch filters 
derived from conventional filters. 

The transient response of notch filters can be kept arbitrarily low 
at large times, such times being defined after the first zero of the 
envelope, t . The time t„ can be kept small by the choice of the number 
of sections m, and/or by choosing /3/<*i large. However, for the interval 
^ t ^ t these methods are not effective. In fact, it has been shown 
above from the initial value theorem (27) that, at t = 0, the envelope 
of the response is unity independent of the filter parameters. A low- 
pass filter combined with a notch filter would cause the transient 
response to be zero at t = 0, and behave, for small t, as t k ~ l if k is the 
order with which the transfer function goes to zero as s — » °° . However, 
with a given low-pass filter the transient response may not be reduced 
to a desired level at small t. An additional method of reducing the 
response by the use of phasing sections is discussed subsequently. 



4.4 Numerical computations 

To illustrate some of the properties of notch filters, numerical 
computations for a three-section filter (?n = 2), with the parameters 
given in Table I, have been performed. Figure 4 shows the filter 
response to a step function and Fig. 5 the response to a stepped cosine 
function. The computed results are essentially in agreement with 
those obtained based on the approximate method. Figure 6 shows the 
transient response of this filter followed by a C-message weighting 
filter, when excited with a stepped cosine function. A comparison of 
Figs. 5 and 6 shows that the C-message filter reduced, as expected, 
the first lobe of the response, affected only slightly the second lobe 
and increased the subsequent lobes. Increasing /3 and hence the 3-dB 
bandwidth of the notch is not very effective in reducing the second 
lobe. This led to the investigation of phasing sections as a means of 
reducing the transient response. 
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Fig. 4 — Unit step response of notch filter. 
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Fig. 5 — Stepped cosine response of notch filter. 
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2.0 2.5 3.0 

TIME IN ms 
Fig. 6— Response of notch and C-message weighting filters to stepped cosine 
function. 

V. TRANSIENT RESPONSE OF PHASING SECTIONS 

The transient response of a phasing section is considered when 
excited by a damped sine function representing the output of the 
notched low-pass filter combination. Let the transfer function of the 
phasing section, T p (s), be given by 

s 2 - cs + d? 



r p( s ) - s2 _j_ cs + d 2 
and the Laplace transform of the damped sine function, F d (s), by 



(28) 



F d (s) = 



s 2 + as + b 2 



^29) 



with Wd = Vb 2 - (a/2) 2 . The Laplace transform of the response, F(s), is 

Wd (s 2 - cs + d 2 ) 



F ^ " s 2 + as + b 2 s 2 + cs + d 2 
Equation (30) can also be written 



(30) 



F(s) = 



U>d 



s 2 + as + b 2 
— 2co)d 



-n — : ; rs ~i _o i i m 



(31) 

s 2 + as + 6 2 ' s 2 + cs + d 2 

where the constants A, B, C, and D are obtained by comparing (30) 
and (31). 
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It can be readily shown that the time response corresponding to 
(31), /(*), is 

f(t) = e-< ffl '«' |~sin u d t - ^ sin (<a d t + O )1 



+ 2cd aa e _ (cmt gin M + 0i) (32) 
r o>i 



where 



and 



r = V(d 2 - 6 2 ) 2 + (c - a)[6 2 c - d 2 a], (33) 

x fl (d 2 - b 2 )u d , OA , 

** «• = Hc _ a) _ a/2(d s _ y) . (34) 

*"*- ,,(.- If r e 7 2 £-*r (35) 



= vRI7 



«i = i/* " o • (36) 



The first term in (32) is the damped sine function and the other two 
terms are introduced by the phasing section. In order that the last 
two terms be of significance, it is necessary that these terms be com- 
parable to the first term. This will be the case for d = b, for which 
(32) reduces to 

f(t) = r "'«' 2jL2 s i n Udt ?c_ ui g _ (c/2)t sin tf (3?) 

a — c a — c o)i 

If, in addition, (a/2) 2 « b 2 and (c/2) 2 « d 2 , (37) simplifies further 
and, for (a/2 — c/2)t <«C 1, (37) is approximately given by 

f(t) « e- a ' 2t sin W(l - c«). (38) 

Equation (38) contains the damped sine function but modified by 
the term (1 — ct). This term can be used to introduce a zero in the 
time vicinity where the damped sine function assumes a maximum 
value. 

To illustrate the above, a phasing section is introduced to modify 
the impulse response of a C-message weighting filter. The computed 
impulse response of the filter is shown in Fig. 7 and has relatively 
large values in the time vicinity of 0.4 ms. The impulse response 
modified by a phasing section is shown in Fig. 8. The phasing section 
parameters are c = 2-10 3 and d = 10 8 . These parameters have been 
chosen on the basis of the above analysis. It is evident the large values 
of the response have been reduced, but the modified response has 
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4.6 5.0 



Fig. 7— Impulse response of C-message weighting filter. 

appreciable values for a much longer time duration than the initial 
response. This behavior can be explained on the basis of Parseval's 
theorem, 14 since the absolute value of the Fourier transform of the 
response is the same with and without the phasing section. 




Fig. 8— Impulse response of C-message weighting filter with phasing section. 
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Table II — Two-section filter with phasing section 
(3-dB bandwidth 390 Hz, 30-dB bandwidth 80 Hz) 





Parameters in kHz Units 




Vnt 


2 


Odi 


2 
Udi 


1 

2 
3 


0.09 
0.09 
-1.8 


40.27 

40.27 

100.00 


1.35 
1.94 
1.8 


27.91 

51.10 

100.00 



Phasing sections can be used to reduce the overshoot of the response 
of a notch filter followed by a low or bandpass filter and excited with 
a stepped trigonometric function at the notch frequency. Such sections 
are of particular importance where the transient response has to be 
modified without affecting the amplitude of the frequency response or 
where a modification is needed at times shortly after the beginning of 
the response. However, such sections may also introduce considerable 
distortions of the unit step response. 

As an example, the performance of a 1010-Hz notch filter with 
and without a phasing section is considered. The transfer function, 
T(z), of the filter and phasing section can be written as 



rw = ft :;+^+i' - 

,=1 Z l + (TdiZ + Udi 



(39) 



The parameters in (39) are listed in Table II. 

This filter was derived from a Tschebyscheff low-pass filter, and 
an operational amplifier version was synthesized and built. Figures 9a 
and 9b show the computed response to a stepped cosine of the filter 
combined with a C-message weighting filter without and with the 
phasing section. Figures 9c and 9d show photos of the corresponding 
oscilloscope displays obtained with the actual filters. Good agreement 
was obtained between the computed and measured response. The 
effect of the phasing section on the response is evident in this figure. 
About a 4-dB reduction in the overshoot was obtained with the 
phasing section. 

VI. CONCLUSIONS 

The transient response of a class of notch filters which are derived 
from low-pass filters by a frequency transformation was investigated. 
General expressions for the transient response due to a unit step 
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nnw 



/A* A,. 



Pi „ Q_Comp»ited and measured response of notch filter with C-message weighting 
filter to a stepped cosine: (a) and (c) without phasing section, (b) and (d) witn 
phasing section. 

function and a stepped trigonometric function have been obtained in 
terms of the low-pass impulse response. 

The transient response of certain types of notch filters can be formu- 
lated approximately in terms of Laguerre functions. These filters have 
been examined in detail and some general properties of the notch 
filter response have been deduced from this formulation. 

Notch filters may considerably distort short time pulses (time 
duration less than one-half notch frequency period). The amount of 
distortion depends on the notch depth. 

The response of notch filters to stepped trigonometric functions can 
be kept at low levels only after a certain time interval from the begin- 
ning of the response. The length of the time interval depends on the 

filter parameters. 

A method of reducing the transient response at short time intervals 
by the use of phasing sections was presented. This method may prove 
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particularly useful in applications where it is necessary to modify the 
transient response without affecting the frequency response. 
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APPENDIX A 
Transient Integrals 

A.I Unit step response 

The relationship between the step response of the notch filter and 
the impulse response of the low-pass filter from which the notch filter 
is derived is given by (6) of the text, 

hs(f) = ~ f v f" ex P (~ ^T7? u ) ^ L ^ du dz - (40) 

After interchanging the order of the integrations and expanding the 
exponential function, (40) can be written 

Equation (41) contains a sum of inverse Laplace transforms of a 

tabulated form 1516 

= T ( l -^ \j2vl2o> ^ut- v*]g(u)du, (42) 
Jo \ o) u / 

where g(u) is the inverse Laplace transform of G(z). 

Equation (42) is of its own interest, since it may be used to obtain 
the step response of a bandpass filter derived from a low-pass filter by 
a frequency transformation. A direct derivation of (42) follows. 

Using the definition of G[z + (ul)/z}, (42) can be written 

- s? L£* /." exp H* + 1) u ] giu)du - (43) 
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After interchanging the order of integration and expanding the 
exponential function in a power series, (43) can be written 

2irj J r z 2v+i \ */ 

= i fW) f ft*- 1 £ ( 7 < ^Ti <tedtt - ( 44 ) 

2irj 7o Jy r^om\z n+2v+l 

The inverse of each term in (44) is zero for a negative exponential 
argument. For a positive argument u ^ t, the inverse is readily ob- 
tained; hence, 

~h g{u) £o iU (m+2,)!^ (45) 

The summation of the terms in (45) gives a Bessel function 17 of order 
2v, and hence (42). Using (42), (41) can be written 

/■oo /•« rt oo ^ /3m) n x n ~ l 

Mt) = / o /*«*.+/, AM /, „?, SrT^ (^Tf! 

• J (2o)„Vxe - x 2 )dz dw. (46) 
The series in (46) can be summed yielding a Bessel function of 
order one ; hence, 

Mt) = r u(u)du - f r f L (u) J^ji{24&v)du 

Jo Jo Jo \ * 

■J (2u o ^tx - x 2 )dx. (47) 

The integration with respect to u can be interpreted as a Hankel 
transform of the low-pass impulse response. For an impulse response 
which is not very oscillatory, and for « «„, the Hankel transform 
will be a slowly varying function in comparison to the Bessel function. 
Under these conditions an approximation to (47) can be obtained by 
using the method of stationary phase. 

A.2 The stationary phase approximation 

The stationary phase method 18 approximates integrals, /, of the 
following type : 

/ = f g(x)eW>dx, (48) 

where A; is large and g(x) is a slowly varying function. The approxi- 
mation considers only contributions from the vicinity of stationary 
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points where (d\f/)/dx = 0, and is of order (1/A;). The approximate 
value of the integral (48) is 

7 *?«W>J^^<->. (49) 

To bring (47) to a form suitable for evaluation with the stationary 
phase method, the Bessel function is expressed in terms of modulus 
and phase. 19 

J {z) = M (z) cos e o (z) (50) 

with 

d (z) =z~\ + ••(«), (51) 

where 5 (z) and M (z) are slowly varying functions for large z with 
lim z -.oo 8 (z) = and lim z _oo M (z) = ^2/(tz). 

With 2 = 2a> Vto — x 2 , the integral in (47) has a stationary point at 
x = t/2. The approximate value of (47), obtained by using (49), is 

fira(f) ^| o °° h(u)dt - y[^ o Mo(<*at) cos U(u> o - |1 

• r hiu) J^ Jii^mWu. (52) 

For large values of t such that M (z) and 6 {t) can be approximated 
with their asymptotic values, 

fsa (0 » r -M«)* - ^^ r fM \l^ Ji(-&M)du. (53) 

JO O)o Jo Y t 

It is of interest to note that the stationary phase method gives the 
correct value for the integral 

f ' J (2co Vu* - w 2 )dw = ^-^ • (54) 

./o wo 

This integral can be evaluated exactly by using (42) with v = and 
g(u) = 1.0. The left-hand side of (42) is readily inverted yielding (54). 

APPENDIX B 

Comparison of Exact and Approximate Solutions 

Consider a notch filter transfer function 



™=(dwhr- (M) 
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The Laplace transform of the time response due to a unit step 
function is 

*«<» -^-ydh^r- (56) 

For m = 0, the time response is 

Us(t) = (l " £e-«" 4)l sin ««)-l(0, (57) 

where 

- - ^-(f)' (58) 

A comparison of (57) with (19) shows that both are of the same 
form but o> is replaced by w . Hence, the approximation is of order 

(/3/"o) 2 - 
For m = 1, using tables of Laplace transforms, 

+ (£Y|8to-W»«C0B«rf. ( 59 ) 

Neglecting terms of order (/3/o>) 2 against one, (59) can be written 
fasit) = 1 - - e-c»'"'Ll (| t) smuot 

+ ( Jt-Ypte-W'eo&aJ, (60) 
\2w / 

where from (16) 

L}(x) = 2 - x. (61) 

The terms multiplied by sin w t are up to order (/3/w ) 2 the same 
as in (19). 

For m = 2, 

«»-!•--' [»+5S+l(S)' 

-?(»+£)+ifi)'( 1 -S)]— 

+ S-"- *[i + i5-?]— '• (62) 
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Neglecting terms of order (/3/oi ) 2 and higher against one yields 
f NS (t) - 1 - L r-C»»i£l ( 1 1) sin t*J 



?2 



.■#r7- 



+ - 2 «-«"» ' ? [7 - 00 cos tfa f, (63) 



u 



8 



where from (16) 



I4(«) = 3- 3x + ^- (64) 



REFERENCES 

1. D. L. Favin, "The 6A Impulse Counter," Bell Laboratories Record, March 

1963, pp. 100-102. 

2. Members of the Technical Staff, Bell Laboratories, Transmission Systems for 

Communications, Fourth Edition, 1970, pp. 32-34. 

3. M.E. Van Valkenburg, Introduction to Modern Network Synthesis, New York: 

John Wiley, I960, pp. 479-486. 

4. P. M. Morse and H. Feshbach, Methods of Theoretical Physics, Part I, New York: 

McGraw-Hill, 1953, p. 944. 

5. W. T. Howell, "On a Class of Function Which are Self-Reciprocal in the Hankel 

Transform," Phil. Mag., 26, Series 7, April 1938, pp. 622-628. 

6. A. Erdelyi, et al., Tables of Integral Transforms, Vol. 2, New York: McGraw-Hill, 

1954, p. 43. 

7. E. A. Guillemin, Synthesis of Passive Networks, New York: John Wiley, 1957, 

p. 489. 

8. Ref. 3, p. 343. 

9. Ref. 6, Vol. 1, p. 175. 

10. I. S. Gradshteyn and I. M. Ryzhik, Tables of Integral Series and Products, New 

York : Academic Press, 1965, p. 1037. 

11. Ref. 10, p. 718. 

12. K. W. Henderson and W. H. Kautz, "Transient Response of Conventional 

Filters," IRE Trans. Circuit Theory, December 1958, pp. 333-347. 

13. M. F. Gardner and J. L. Barnes, Transients in Linear Systems, New York: 

John Wiley, 1951, p. 267. 

14. Ref. 4, p. 456. 

15. Ref. 6, Vol. 1, p. 133. 

16. N. W. McLachlan, Modern Operational Calculus, London : The MacMillan Com- 

pany, 1949, p. 45. 

17. Ref. 10, p. 951. 

18. A. Erdelyi, Asymptotic Expansions, New York: Dover Publications, 1950, pp. 

51-57. 

19. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, New 

York : Dover Publications, 1965, p. 365. 



304 THE BELL SYSTEM TECHNICAL JOURNAL FEBRUARY 1974 



